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Abstract 


A  coherent  system  with  independent  non-renewable  components  with  arbitrary  lifetime  dis¬ 
tributions  is  considered.  A  simple  observation  leads  to  a  hierarchy  of  upper  and  lower  bounds 
that  converge  to  the  exact  system  reliability.  The  simplest  of  these  bounds  is  shown  to  be 
tighter  than  the  bounds  of  Gertsbakh  (1985). 
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Abstract 


1.  Introduction 


Consider  a  coherent  system  consisting  of  n  independent  non-renewable  components  {1,  2,..,  n], 
with  lifetime  distributions  »-  1,  2,..,n.  Let  Jf-  {*(/),'>  0}  be  the  vector  performance 
process  of  these  components  such  that  X[t)  ®  (Xi(t)tX2(t),..tXn(t))  and  A$(/)  takes  the  value 
0  if  component  i  is  up  at  time  t  and  1  otherwise.  The  state  space  S  of  X  is  then  {0,  1}". 
Let  $  be  the  coherent  structure  function  of  the  system.  That  is  $:S  -»  {0,  1},  is  non-decreasing 
on  S  and  when  <f>(X(fl)  is  0  the  system  functions  at  time  t  and  is  1  otherwise.  The  reliability 
R(t)  of  the  system  at  time  t  is  then  given  by  1  -  £{£(*(/))}.  It  is  of  interest  to  obtain  R. 
In  general  it  is  not  possible  to  obtain  explicit  formula  for  R  and  one  resort  to  algorithmic 
methods  to  compute  R.  Except  in  some  special  cases  [e.g.  Agrawal  and  Satyanarayana  (1984), 
Provan  and  Ball  (1984),  and  Shanthikumar  (1982,  1984)],  the  problem  is  in  general  NP-hard 
[e.g.  Buzacott  (1980),  and  Ball  (1984)].  Consequently  approximation  and  bounds  have  been 
developed  [e.g.  Ball  and  Provan  (1983),  Gertsbakh  (1985),  Shanthikumar  (1984),  and  Shogan 
(1976)].  In  this  paper  we  derive  alternative  bounds  for  R  and  compare  them  to  that  of 
Gertsbakh  (1985). 

Let  h(p)  a  £{$Uf(0)}.  where  p,  a  £{AJ(0}  -  £,(/),  and  let  ((?,£)  be  a  partition  of  S  such  that 
Gmm  -0,xe  S}  and  define  £//  — (x:  |  |x|  | S},/  — 1,2, where 

|  |  x  |  |  a  X  x,.  Suppose  min  {  |  |  x  I  |  :x  €  B]  —  r.  Then  where 

VjwmB  fl  l/-'j  —  r,  r  +  1, ...  n.  With  these  observations  one  has,  *“* 

n4|n<i-'’4-  (1) 

-  .  (2) 

t-r 

where 

^-s{n4{n<>-4  •  <« 

WXi)  “  {kutfc  —  O.k  — 1,2,  ..,n},  is  the  set  of  working  components  in  state  x,  and 
iP(x)-  {l,2,..,/i}- W(2). 

Suppose  we  have  upper  and  lower  bounds  Bf  >  A/  >  Cf,  f  —  r,r  +  1  ,..,n.  Then  we  will  have 
a  sequence  of  upper  and  lower  bounds  for  h(p)  and  consequently  for  R.  Specifically 


Bounds  for  System  Reliability  Function 


1 


(4) 


/-*■  /-«.  “  f— r 

» 

for  m=r,r+lt..n+l,  [  we  use  2y,~0,  h<a  far  all  yj.  Note  that  the  above  bounds  get  tighter 
as  m  increases  and  give  the  exact  unrealibflity  when  m  *  n  +  1.  This  allows  one  to 
progressively  compute  Ar,Ar+\, ....  until  a  desired  accuracy  is  found. 

In  Section  2  we  develop  some  upper  and  lower  bounds  for  At  and  compare  it  to  that  of 
Gertsbakh  (198S).  A  simple  numerical  example  is  given  in  Section  3. 


2.  The  Bounds 

Without  loss  of  generality  assume  that  the  components  are  numbered  so  that  p\  >  p2  £  ...  £  p„- 
Since  p;  *■  /}(/)  depend  on  t,  this  numbering  may  have  to  be  altered  for  different  values  of  t, 
unless  /=i(/)  >  F2(t)  >  ...  £  fn(0,  for  all  /  £  0. 

Proposition  2.1 :  Let  |  |  Vf  |  |  be  the  number  of  elements  in  the  set  ty.  Then 

1 1  r>  1 1  n  a = c}« 

(.  i-l  J  IK+I  ) 

jn*J[na  -vj  ® 

Proof-.  It  is  easily  seen  that  for  any  xc  Vf, 

n <>  s  { n  4{  n <> -^} 

l. J  j  'ttir(r)  *  v«l»Tx)  * 

‘{B-EdH 

The  result  now  directly  follows  from  (3). 
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To  use  the  bounds  and  one  needs  to  know  |  \  Vg\  I .  Once  |  \V(\  | —  r,  r  +  1, .., 
are  known  one  may  use  the  recursion  -{11  V/+i  I  1  /  1  1 V/ 1  jj  \pt+ 1  /  (1  - Pt+i)}  to 
compute  B^\  t  —  r,  r  +  1, ...  [A  similar  recursion  may  be  naed  for  Q  *  J  —  r,  r  +  1, ...]  In  some 
cases  obtaining  1  |  V/  |  |  itself  may  turn  out  to  be  computationally  difficult.  In  such  a  case 
the  observation  that  I  I  ^  I  I  <  f  ^  j  leads  to  a  modification  of  2?^  as  follows: 


Proposition  2.2 : 


»+i-< 

We  will  next  obtain  a  modification  of  Bj2\  Define  Pf"  —  2  pj/(n+l  —  i)  and 
Pi  —  2 pi/  (n  +  1  —  0,  i  —  1, 2, ....  n.  For  any  y  e  S  let  >"1 


2,00 


Pt  ^ 

\-p-  if  yf  —  O' 


Then 


Proposition  2.3 :  For  any  y  e  U/, 


(7) 


jP/oo/:  Let  P  be  the  collection  of  all  permutations  of  {1,2,..,  n}.  Then 


(8) 


where 
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Now  consider  a  subset  P n**1  of  P,  such  that  the  first  (n-1)  elements  in  each  of  the  permutations 
in  Pn~l  are  identical.  Then  for  any  v'  c  Pn~i 

2  2 

wtyt~l  *-l  <-l  9tr"~l 


l-l  rtr*-1 

-  2 

since  by  definition  Qniy)  >  4,(w,y)  for  all  v  e  P"-1.  Note  that  |  I/*"-1 1  I  -  1.  Therefore 

2  fl*. (*.>-) s  2 

**r  f-i  wp  i-i 

Now  defining  Pn~ 2  to  be  a  subset  of  P  such  that  the  first  n-2  elements  of  every  permutation 
in  there  is  identical  one  sees  that 

m-2  n—2 

mtf»-2  (-1  wtr”~2  1-1 

since  by  definition  2. Q„_i Cy)  >  2  gn-\{v,y).  Note  that  I  I  Pn~2  11-2.  Then 

•cl*-2 

2  ^  2  Tbw  fl 

9€P  1-1  9€p  1—1  J~n—  1 

Continuing  this  way  and  observing  that  |  |  P  |  |  -  n\.  Equation  (8)  leads  to  the  required  result. 

□ 
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Since  (7)  may  give  different  upperbounds  for  different  / s,  it  is  appropriate  that  we  choose 
the  mim'mum  of  such  bounds.  In  this  regard  define  v  *  to  be  a  permutation  of  {1,2 
such  that 


p+ 

^•(l) 


PZ» 


(2) 


pZ, 


1  -  ^V(l)  1  ~  ^W*(2) 


:r—  £  •—  < 


(n) 


1  -  P~. 


(9) 


(«) 


Then  using  a  straightforward  pairwise  interchange  it  can  be  shown  that 
Proposition  2.4 


a{oc,w} 


rp*«  fi (i  -  w- 

/-I  /_/+! 


Equivalently 


'  *  (-i  j-t+i 


(10) 


l-k 


We  will  next  provide  an  alternative  upper  bound  for  2^4/  which  is  tighter  than  that  given 
by  Gertsbakh  (1985). 

Proposition  2.5: 


t-k  N  /  f-1  f-k 

Proof:  Observe  that 

n 

£P{at  least  k  components  are  down  at  time  t}. 

t-k 

Since  for  any  ir  e  P, 


(11) 


(12) 
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P{components  ir(l),ir(2),..,ir(fc)  are 


it 


down}  -  iw 

i-i 


one  has  from  (12) 


^'Sp^hjrSri'V.v  <13> 

t-k  K-'-n  K>’  wtf  i-i 

Using  an  analysis  similar  to  that  used  for  the  proof  of  Proposition  2.3  one  gets 

2  rKo  fK  <14> 


w  €  P  l-l 


Combining  (13)  and  (14)  one  obtains  the  desired  result. 

□ 

When  the  lifetime  of  component  i  is  exponentially  distributed  with  rate  A,-,  one  has 
Pi  -  1  -  exp(— A,/),  i - l,2,.../i.  With  Aj  £  A2  £ ...  £  A„, 


■+i— i 


: i — 

n  +  l-i  11 


where  A,  =  X  Ay.  Therefore  from  (11),  one  gets 


n  .  .  r+l 


(15) 
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The  right  hand  side  of  (15)  is  the  bound  obtained  by  Gertsbakh  (1985).  Clearly  (11)  is  a 
tighter  bound  and  is  applicable  to  any  component  lifetime  distributions. 

In  summary,  one  can  use  (4)  along  with  (5),  (6),  (10),  and  (11)  to  obtain  alternative  bounds 
for  the  system  reliability  function  R.  Since  the  computation  of  these  bounds  are  relatively 
simple,  all  these  bounds  may  be  computed  and  the  best  value  can  be  used. 


3.  A  Numerical  Example 


Consider  the  network  shown  below  in  Figure  1.  Suppose  we  are  interested  in  the  reliability 
that  node  s  is  connected  to  node  t  by  at  least  one  path  of  working  edges. 


Figure  1:  Directed  Reliability  Network 


i 

1 

2 

3 

4 

5 

.1 

.1 

.08 

.08 

.08 

Table  1:  Edge  unreliabilities  (/>,). 

The  edge  failure  probabilities  /?,-  =  1,2, 3, 4, 5  are  given  in  Table  1.  Then  the  minimum  number 
of  failed  edges  needed  for  a  cut-set  is  r=2,  and 

*4-  {(0,0, 1,1,1),  (1,1, 1,0,0),  (0,1, 1,1,0)} 

V3  -  {(0,0,0, 1, 1),  (0,0, 1,0, 1), (0,0, 1, 1,0), 
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(0,1, 1,0,0),  (1,0, 1,0,0),  (1,1, 0,0,0), 

(0,1,0, 1,0),  (1,0, 0,0, 1,1)} 

*4=  {(0,0, 0,0,1),  (0,0, 0,1,0),  (0,0, 1,0,0),  (0,1, 0,0,0),  (1,0, 0,0,0)} 


V5-{  (0,0, 0,0,0)}. 

with  |  |  V}  |  |  **  3,  |  ]  1^3  |  |  =»  8,  |  |  Kj.  |  |  =5  and  I  I  Vs  |  |  -» 1.  Exact  values,  upper  and  lower 
bounds  of  A(,l  =*  2,3, 4,5  are  calculated  using  Equations  (3),  (5),  (6),  (10),  and  (11)  and 
tabulated  in  Table  2.  In  this  example  Pf  is  0.088,  0.09,  0.09333,  0.1,  0.1  and  P~  is  0.088, 
0.085,  0.08,  0.08,  0.08  for  i  =  1,  2,  3,  4  and  5  respectively.  Hence  v  *  =»  <  1,2, 3,4,5  >,  thus 
resulting  in  the  same  values  for  the  bounds  and  [Equations  (6)  and  (10)].  Substituting 
the  upper  and  lower  bounds  for  Ai  in  Equation  (4)  one  obtains  the  bounds  for  the  unreliability 
h{p)  [shown  in  Table  3]. 
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m 

2 

3 

4 

5 

£<7(3)  £  >1/ 

f — m 

0.02348992 

0.00483968 

0.00027392 

0.00000512 

£,7(5)  £  Cj!> 

r—m 

0.01786112 

0.00355328 

0.00023552 

0.00000512 

£,7(5)  £ 

t—m 

0.02907712 

0.00571648 

0.00029952 

0.00000512 

Eq{ 6)  £  £P 

f— m 

0.08493952 

0.00707072 

0.00029952 

0.00000512 

£<7(10)  £ 

/— m 

0.08493952 

0.00707072 

0.00029952 

0.00000512 

Eq(  11)  £  £*4) 

/— m 

0.0792 

0.007392 

0.0003696 

0.000007392 

Table  2:  Exact  Values,  Upper  and  Lower  Bounds  for 

f— m 


m 

Lower  bound 

using  Eq(5) 

Upper  bounds  Using 

Eq(5) 

Eq(6) 

Eq(10) 

Eq(ll) 

E 

0.01786112 

0.02907712 

0.08493952 

0.08493952 

0.0792 

D 

0.02220352 

0.02436672 

0.02572096 

0.02572096 

0.02604224 

B 

0.02345152 

0.02351552 

0.02351552 

0.02351552 

0.0235856 

1 

0.02348992* 

0.02348992* 

•0.02348992* 

0.02348992* 

0.0234901472 
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